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Abstract 

We introduce a new regularization of the rotational Navier-Stokes equations that we call the 
Rotational Approximate Deconvolution Model (RADM). We generalize the deconvolution type 
model, studied by Berselli and Lewandowski [5], to the RADM model with fractional regulariza- 
tion where the convergence of the solution is studied with weaker conditions on the parameter 
regularization. 
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1 Introduction 

Let us consider the rotational Navier-Stokes equations in a three dimensional torus T3, 

V-v = 0, (1.1) 
v^t-v xW XV - i^Av + VP = f, (1.2) 

with initial data v{x,0) = vq{x). Here, v is the fluid velocity field, P is the Bernoulli scalar or 
the dynamic pressure P = p + ^u^, / is the external body forces and u represents the viscosity. 

When the flow is turbulent it is well known that equations p.l|) - (|1.2|) are unstable in nu- 
merical simulations. Therefore, the numerical turbulent models are needed for real simulations 
of turbulent flows. The alpha models is an example among these numerical models that have 
attracted a lot of interest from both the pure and applied mathematicians. These models are 
based on a filtering or averaging obtained through the application of the inverse of the Helmholtz 
operator 

A:=I-a'^A. (1.3) 

Similar to the alpha models [ZKIIKiaElEKISlliaiSllES], Olson and Titi ^ used a more 
general filtering with fractional regularization given by 

V = Aqv := (l + a^^i-Af) v for some fixed 6* > 0, (1.4) 



such that V denotes the unfiltered velocity and v denotes the filtered one. This filter is a 
differential filter [12j, that commutes with differentiation under periodic boundary conditions. 

The Approximate Deconvolution Model of turbulence (ADM) belongs also to the class of 
numerical turbulent models. A complete study of the ADM models has been carried out during 
the last years thanks to their well-posedness and their good performance in numerical simulation 
[U \20\ 121] . The ADM model consists in replacing the nonlinear term in the Navier-Stokes 
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equations V • (i? (X" i?) by V • {D^v ® D^v), where the exact Helmholtz operator is used 
to define the deconvolution operator Djsj of order N as fohows: 

N 

Dr, = Y.{I-A-y. (1.5) 

1=0 

Existence, regularity and uniqueness of a solution to this model for a general deconvolution 
of order A^, were proved in [U]. The case = was studied in details in [TU [TS]. Recently, 
Berselli and Lewandowski [5] studied a generalized ADM model that consists in replacing the 
deconvolution operator by a more general deconvolution operator D^^q, with fractional 
regularization. This deconvolution operator is defined by 

N 

DN,e = J2(^-A,'y. (1.6) 

i=0 

For > |, Berselli and Lewandowski ^ showed the convergence of the ADM solution, with 
fractional regularization, to a solution of the average Navier-Stokes Equations when N goes to 
infinity. The existence and the uniqueness of a solution to this model in the case A^ = and 
> ^ are studied in [3]. 

In this paper, we introduce a new model called the Rotational Approximate Deconvolution 
model of turbulence (RADM) . Let v be the filtered or the averaged velocity field which is given 
by the relation (|1.4|) . In order to obtain the RADM, we first write the Reynolds averaged 
rotational Navier-Stokes equations 

v^t - V x\/ X V - i^Av + VP = 7, (1.7) 

which are not in closed form due to the averaged nonlinear term v x V x v that we cannot 
write in terms of v alone. 

As the main goal of turbulence modeling is to derive simplified and computationally real- 
izable closure models, we use the deconvolution operator D^ q to approximate the velocity v 
with Dn^qv and write the equation for v as bellow 



v,t - Dn,9V X V X DMfiV - lyAv + VP = /. (1.8) 
This yields the initial value problem 

V-w = 0, (1.9) 

w^t — D]\f^0W X V X Dn£W — V/S.W + Vq = f, (1-10) 



with initial data w{x, 0) = wq{x) = vq, considered in (0, T) x T3 and subject to periodic bound- 
ary conditions with mean value equal to zero. Note that the couple {w, q) is the approximation 
of {v,P). 

When A^ = 0, the operator D^ g = I. Thus the zeroth order RADM is given by the following 
equations 

V-w = 0, (1.11) 
w^t — w X \7 X w — uAw + \7q = f. (1-12) 
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The model p.9p - (|l.lU|) is a deconvolution model for the rotational form of the nonlinearity. 
Therefore calling it the RADM model seems natural. This model is interesting as it can be 
discretized using methods similar to those used for the rotational form of the Navier-Stokes 
equations. Recent work on the discretizations of the Navier-Stokes equations [1] revealed that 
solving the NSE in this form can lead to better efficiency, especially in the case of low viscosity. 

The purpose of this paper is to study the RADM model with fractional regularization 6 
and look for the limiting case where we can prove global existence and uniqueness of regular 
solutions {wN,qi\f) for a fixed A^. The analysis of the zeroth order RADM ()l.lip - ()1.12p is the 
same as the general case. For this reason, the regularization parameter 9 does not depends on 
A^. Another goal of this paper is to study the question of the convergence of the RADM solution 
when N goes to infinity. Namely, we show that when N tends to infinity and | < ^ < 1 the 
unique solution (wN,qN) of the RADM model satisfies 

V-w = 0, (1.13) 
w^t — Agw X V X Agw — uAw + Vq = f, (1-14) 
w{x,0) = wo{x) = vo. (1-15) 

This result of consistency and stability illustrates how the system approximates the rota- 
tional Navier-Stokes equations. 

The rest of the paper is organized as follows. In section 2, we start by giving some mathe- 
matical tools such as the space of functions we are using and the deconvolution operator D^^q. 
Section 3 provides existence and uniqueness results in the case where < A^ < oo. The last 
section provides the main result of this paper which concerns the study of the convergence of 
the RADM solution when N goes to infinity. Note that the results of this paper require weaker 
conditions on the parameter than that required in [5]. 

2 Mathematical setting 

2.1 Notations 

In this section we fix notation of function spaces that we shall employ. 

We denote by U'iT'i) and H^iT'i), s>— 1, l<p<oo, the usual Lebesgue and Sobolev 
spaces over T3, and define the Bochner spaces C(0, T; A), L^(0, T; X) in the standard way. The 
Sobolev spaces H" = H^{T?j)^, of mean-free functions are classically characterized in terms of 
the Fourier series 

H'=\v{x)= Cfce*-,(cfe)* = c_fc,Co = 0,||^||2 2= ^ \k\^'\cu\^ < ^ 
y keTz keTa 

where (c^)* denote the complex conjugate . In addition we introduce 

H';^ = {veH'; V-^; = Oin T3} , 
H ' = [H')' , = Ll = Hi. 

Throughout we will use C to denote an arbitrary constant which may change from line to line. 

2.2 The generalized deconvolution operator 

The deconvolution operator considered in this paper is a generalized deconvolution operator 
that is constructed by using the Helmholtz equation with fractional regularization [5, 2]. Let 
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a > 0, s > —1, < < 1, V €^ H'^ and let v G fj^+'^s \)q the unique solution to the equations 
(compare to [T7] ) 

a'^\-Afv + v = v, (2.1) 
V ■v = V ■v = 0. (2.2) 

We also shall denote by Aq the operator 

^s- ^^„2e(_^)'e^ + ^_ (2-3) 
The non-local operator Aq is defined through the Fourier transform 

= (l + a^^lfcp^) v{k). (2.4) 

Therefore, one has 

v = Ag^v, (2.5) 

and 

\\v\U+2e < -^\\v\\s. (2.6) 

Let us consider the operators 

N 



i=0 

When 9=1 and for a fixed > 0, we recover the van Cittert deconvolution operator 
used in Large eddy simulation (LES) by Stolz and Adams [T] see also [S]. This deconvolution 
operator, for different values of 9, was used in |10j to study the rate of convergence of the ADM 
model to the mean Navier-Stokes Equations. 

A straightforward calculation yields 

(\ / / 29\i,\2e \ 



Thus 



I E = E A^^(fc)cfce^''■^ (2.^ 

,fceX3 / fceXa 



where we have for all A; G X3, 



DoAf') = 1, (2.9) 
1 < Av^(fc) < iV + 1 for each iV > 0, (2.10) 
and Av^(fc) <Ag:=(l + a'^^\k\'^'^) for a fixed a > 0. (2.11) 



One can prove the following Lemma (see in [5] ): 
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Lemma 2.1. For all s > —1, k G X3, and for each N > there exist a constant C > such 
that for all v sufficiently smooth we have 

\\v\\s,2 < \\Dn,9 (v) \\s,2 <{N+ l)\\v\\s,2, (2.12) 

\\v\\s,2 < C\\DN,e (v) \\s,2 < C||AeiL»|^g(^;)||,,2, (2.13) 

||Ae^L»|^g(U)||,,2 < ||^||.,2, (2.14) 

Ms+e,2 < Cia)\\Ae-2Dl/v%,2. (2.15) 

3 The RADM model: existence and uniqueness of a weak solu- 
tion 

Let us consider the following Rotational Approximate Deconvolution Model (RADM) 

V-w=0, (3.1) 

W^t — DjqeW X V X Djy gW — i^Aw + Vq = f, (3.2) 

[ w = 0, I q = 0, (3.3) 

w{x,0) = wo{x) = v^, (3.4) 



with periodic boundary conditions. 

Definition 3.1. Let f € L'^{0,T; H~^) be a divergence free function and Vq G L^. For any 

< 9 < 1 and < N < oo we say that the couple {w,q) is a " regular" weak solution to 
\3.1\ )- [3^^ if the following properties are satisfied: 

w€C{0,T;H',)nL\0,T;Hl+'), (3.5) 

w,t(^L\0,T;H''-l), (3.6) 

q£L^{0,T;H^'-^T3)), (3.7) 



the couple {w,q) fulfill 
rT 



/ {w^t, (p) - {DNfiW X V X D^fiW, (f) + u{Vw, V(f) + {Vq, (f) dt 
Jo 

= [ {f,ip)dt forallipeL\0,T;H^). 
Jo 



(3. 



Moreover, 

w{0) = wq. (3.9) 



Note that the terms {w^t,v), {Dn,9W x V x D]y^gw,(p) and {Vq,ip) in ()3.8p are defined in 
the sense of duality pairing on H^^~^ x H2~^^. The term h'{Vw,V(p) is defined in the sense 
of duality pairing on x H . The last term (/, (p) is defined in the sense of duality pairing 
on H'-^' X H^'-' 

Remark 3.1. The notion of "regular weak solution" is introduced in f^. Here, we use the 
name "regular" for the weak solution since the velocity part of the solution w does not develop 
a finite time singularity. 
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The main result of this section is to prove the existence of a unique "regular" weak solution 
to the RADM model. 

Theorem 3.1. Assume f € L'^{0,T; H~^) be a divergence free function and vq G L^. Let 
> ^ and let < N < oo is given and fixed. Then problem i3.1\) - (3^\ ) has a unique regular 
weak solution. 

Proof of Theorem 13.11 The proof of Theorem 13.11 follows the classical scheme. For sim- 
plicity we restrict ourselves to the critical case = |, and we drop some indices of 9. We will 
consider only the case when < < cxd and the case = is deduced by replacing D^ g 
by I. Thus we will write ^^Di\f" instead of ^^D]\]£^ and "A" instead of "Ag" expecting that no 
confusion will occur. We start by constructing approximated solutions via Galerkin 

method. Then we seek for a priori estimates that are uniform with respect to n. Next, we passe 
to the limit in the equations after having used compactness properties. Finally we show that 
the solution we constructed is unique thanks to Gronwall's lemma |22j . 

Step l(Galerkin approximation). Consider a sequence {t^''}^]^ consisting of L^-orthonormal 
and if ^-orthogonal eigenvectors of the Stokes problem subjected to the space periodic condi- 
tions. We note that this sequence forms a hilbertian basis of L^. 
We set 



N 

w^{t,x) = yc^{t)^^{x), 



r=l 



ik-x 



and q''{t,x) = ^ g^(t)e 



We look for {w"'{t,x),q'^{t,x)) that are determined through the system of equations 

= {f,^n, r = l,2,...,n, 

and 



(3.10) 



(3.11) 



Ag" = V • n" (L'at'u;" X V X DnW) . (3.12) 

Where the projector n"" assign to any Fourier series Qk^^^^ the following series g^.e 

fcez»\{o} fcez3\{o},|fc|<n 
Moreover we require that satisfies the following initial condition 

N 

^"(0,.) = ^o" = E^o¥'^'(^)' (3-13) 

r=l 



ik-x 



and 



1 



Wq — )• Wq strongly in when n — )■ cx). (3-14) 



The classical Caratheodory theory |23) then implies the short-time existence of solutions to 
()3.1ip - ()3.12p . Next we derive estimates on c" that is uniform w.r.t. n. These estimates then 
imply that the solution of (|3.1ip - (j3.12|) constructed on a short time interval [0, T"[ exists for all 

t G [o,r]. 
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Step 2 (A priori estimates) Multiplying the rth equation in (j3.1ip with as |fc| 3L)^c"(t) + 
Dj\fC^{t), summing over r = 1,2, ...,n, integrating over time from to f and using the following 
identities 

{w^t,ADNW^) = il||Ai4^"||2 (3.15) 



1 



1 



(-Ai(;",AD^k;") = \\A^ Dj^w^'Wl^, (3.16) 
(7, AD^w^) = (A5Z)|7, a'^D^w^), (3.17) 

and 

{Dnw^^ X V X DnW^, ADj^w"^) = [DnW^ x V x DnW^, D^w'^) = 0, (3.18) 
leads to the a priori estimate 



t . 1 



i„.i„i „„o r„.i 

2 



/ {A^DlJ,A^Dlw'') ds + \\\A^Dlv^\\l 
Jo ^ 



(3.19) 



Using the duality norm combined with Young inequality and inequality (|2.14p we conclude from 
that 

sup \\A-2Dlw'^\\l + u f \\A^DlwXi,2 ds < \\v^\\l + - T \\f\W2 ds (3.20) 
te[o,T"[ JO ^ Jo 

that immediately implies that the existence time is independent of n and it is possible to take 
We deduce from iK2U\\ and (l2J5]l that 

G L°°(0,r;/i"|) nL2(0,T;J/|). (3.21) 

Thus it follow from j^TT^ that 

Dnw" G L~(0, T;H^)n L^{0, T; Hi), (3.22) 
From ()3.22p and by using Holder inequality combined with Sobolev injection we get 

Dnw"" X V x Dnw"^ G L2(0, T; H^). (3.23) 
From ([3:23|) and dM} it follows that 



L>Arw" X V X Datiu" G L^(0,r;^nr). (3.24) 
Consequently from the elliptic theory eqs ()3.12p implies that 

/ dt < K. (3.25) 

Jo 6'^ 

From eqs. dXTT]) . ([Mil), dST^ and ([H:^ we also obtain that 

T 

\\w'\\\\^dt<K. (3.26) 
' 6' 
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Step 3 (Limit n oo) It follows from the estimates (|3.21|) - (|3.26|) and the Aubin-Lions 
compactness lemma (see [19j for example) that there are a not relabeled subsequence of (if", g") 
and a couple {w,q) such that 



w'' - 


w 


weakly* in L°°{0,T; H^), 


(3.27) 


Dnw"" - 


Dnw 


weakly* in L°°(0, T; 


(3.28) 


w"" - 


w 


weakly in L'^{0,T; nj), 


(3.29) 


Dnw"" - 


Dnw 


weakly in L'^{0,T; hI), 


(3.30) 




w,t 


weakly in L'^{0,T; H^l), 


(3.31) 




Q 


weakly in L'^{0,T; H^T^)), 


(3.32) 


- 


w 


strongly in L'^{0,T; H^), 


(3.33) 


Dnw"" - 


Dnw 


strongly in L^{{),T; Hi), 


(3.34) 



From ([330]) and (|33i]) it follows that 

DnW^ X V X DnW^ Dnw X V X Dnw strongly in L^(0, T; (¥3)^), (3.35) 



Finally, since the sequence Dnw'"- x V x D^w"- is bounded in L {0, T;H e (Ts)'^), it converges 
weakly, up to a subsequence, to some x i^i -^^(0) H~"e (T^)^) . The result above and uniqueness 
of the limit, allows us to claim that x = D^w x V x Dj^w. Consequently 

x V X DnW Dnw x V x D^w weakly in L^{0,T; H^I (Tsf), (3.36) 

The above established convergences are clearly sufficient for taking the limit in (j3.1ip and for 
concluding that the velocity part w satisfy ()3.8p . Moreover, from ()3.29|) and ()3.3ip one we can 
deduce by a classical argument ( see in [2j) that 

w eC{0,T;Hl). (3.37) 

Furthermore, from the strong continuity of w with respect to the time with value in He we 
deduce that w{0) = wq. 

1 1 5 
Let us mention also that D^w + {—A)^ Dj\jw G L^(0, T; ), hence ADj^w is a possible 

1 1 

test function in the weak formulation ()3.8p . Thus A^D^w verifies for all t G [0, T] the following 
equality 

1, 



^\\a'2dI,w\\1 + u f \\A^Dlw\\l^ ds 



/' 

JO 





. 1 i , . 1 



{A^Dlf^A^Dlw) ds+ ^-WA^DlvoWl 



(3.38) 



2' 



Step 4 (Uniqueness) Since the pressure part of the solution is uniquely determined by the 
velocity part it remain to show the uniqueness to the velocity. 

Next, we will show the continuous dependence of the solutions on the initial data and in 
particular the uniqueness. 

Let ^ > I and let (lOi, qi) and (w2, 0.2) be any two solutions of ()3.ip - ()3.4p on the interval [0, T], 
with initial values if i(0) and 1^2(0). Let us denote by 5w = W2 — W1. We subtract the equation 



$LaTeX: 201 3/3/2$ 



for wi from the equation for W2 and test it with ADj^{5w). We get using successively the fact 
that the averaging operator commutes with differentiation under periodic boundary conditions, 
the norm duahty, Young inequahty, Lemma l2.1l 



~\\A-2Dl,6w\\l + u\\VA^Di6w\ 



< {DNiw2) X V X Dm{w2) - Dn{wi) X V X Dn{wi),AD^{6w)) 
< {DnW2 X V X DnW2 — Dnwi X V X DnWi,Dn5w) 



(3.39) 



6'^ 2 6 



1 

2 



< —^-^ —\\DnSw X V X DnWi\\_7 „ + -—\\5w\\i „. 

By using ()2.15p we have 



1 (i 1 2 i 11,-^ r Il2 

-—as \\ow\. 1 + z^as || Voiull 1 

2 at " 



6 



1 d 



1 



1 



(3.40) 



< ~WDl8wf2 + v\VA-^Dlbwf2 



From ()3.39p p.40p and using ()2.12p . Holder inequality and Sobolev embedding theorem we get 

^ ille Il2 lll^c- Il2 

—as llotulll 2 + z^a^ II vowll^ 
, C(a)(jV + l)6 ,„ ,,2„ ' ^^-^^^ 

V 6 6 

Using Gronwall's inequality, since ||ifi||7 € L^([0, T]), we conclude the continuous depen- 

6 

dence of the solutions on the initial data in the L°°([0, T], iJ^) norm. In particular, if bw^ = 
then 5w = and the solutions are unique for all t G [0, T]. Since T > is arbitrary this solution 
may be uniquely extended for all time. 
This finishes the proof of Theorem 13.11 



4 Limit when N ^ oo 

Let {w]\f,q]\f) be the solution of the RADM model constructed above. In this section, we 
take the limit N ^ 00 and look for the optimal value of 9 in order to show that {wi\i,qi\f) 
converges, up to subsequences, to a solution of the mean rotational Navier-Stokes equations. 

Theorem 4.1. Let {wnjQn) be the solution of i3. l\) - (3^\ ) for a fixed a > 0, N > and 
^ < 9 < 1. There is a subsequence Nj —)• 00 as j —> 00 such that: {w]y.,q]y.) — > {w, q) as j —)• 00 

where {w,q) G L'=^{[0,T]; H^) n L'^{[0,T]; H^+^) x L^{[0,T]; H'^+'^^iTs)) is a distributional 
solution of the following system with periodic boundary conditions 

V-w = 0, (4.1) 



w^t — Agw X V X Agw — uAw + Vq = /, (4.2) 



w = 0, q = 0, (4.3) 
w{x,0) = wo{x) = vo. (4.4) 

The sequence w^. converges strongly to w in the space L^([0,T]; H'^) for all s < 1 + 9 , while 
the sequence qat. converges weakly to q in the space L^([0, T]; i7^2+2^(T3)). 
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Remark 4.1. (1) The restriction in the case of the ADM model in is 6 > ^. Our result is 
sharper for the RADM model and it should he clear to the experts that the same approach used 
here can he used to the ADM model to get similar results. 

(2) For 6 > 1, one can uses the same method developed in to show the convergence of RADM 
model to the mean rotational Navier-Stokes equations. 

In order to prove Theorem 14.11 we need to reconstruct a uniform estimates for wisf. Since 
some of the a priori estimates used in the proof of Theorem [3?T] dependent on we can not use 
them. Instead we show in the proof of Theorem 14.11 that Dj^^eWi^ belongs to L°°([0, T]; Z^) n 
L^([0, T]; iJ"^) uniformly with respect to N. 

Before proving Theorem 14.11 we first record the following Lemma. 
Lemma 4.1. Let ^ > and assume thatv e L'^{[0,T],H'^^). Then 

DnmV AqV strongly in L^(0, T; L^), (4.5) 
and there exist a constant C independent from N such that 

\\DN,ev\\2 < W^evh < C{a^)\\v\\2e,2- (4.6) 

Proof. The first part of this lemma is given page 18 in [5]. The second part is a direct 
consequence from the propriety (|2.1ip of the operator Dj\[^q, the fact that 

and Poincare inequality. 

Proof of Theorem 14.11 The proof of Theorem 14.11 follows the lines of the proof of the 
Theorem 4.1 in that we have to modify in order to treat the cases when < |. Thus we 
will use [5] as a reference and only point out the differences between their proof of convergence 
to the mean Navier-Stokes equations and the proof of convergence in the present study. First, 
we need to find estimates that are independent from N . 



Step 1 (Uniform estimates with respect to N) Following the proof of Theorem 13. II we obtain 
that the solution of (I3.1I1-(I3.4I) satisfies 



1 i 1 /■* i i 



1 . 1 (4-7) 

^" • 2 7-12 :77„||2 



Thus 



11 r 1 1 If 

sup \\AlDjfgWN\\l + iy \\A^D^^gWN\\l2 ds <\\vq\\1 + - \\f\\\2 ds (4.8) 
ie[o,T] ' Jo ' ^ Jo 

consequently we can bound the right hand side by a constant C which is independent from N . 
We deduce from (14.81) and Lemma |2. II that 



DNfiWN G L°°(0,r;l.^) n L'^{0,T;HI), uniformly with respect to N, (4.9) 
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and 

WN e L°°(0, T; i?^) n l2(0, T; H^"^'^), uniformly with respect to N. (4.10) 
We observe from (|4.9|) that 

DnA^n) X V X Dn^wn) e L\0,T;H-'2{T3f). (4.11) 
Thus from ()4.1ip and (|2.ip we obtain 

DnA-^n) X V X DnA^n) G L2(0,T;/7-t+2«(T3)=^). (4.12) 
For the pressure term g^r, we deduce that it verifies the fohowing equation 

A^AT = V • {DNfiWN X V X Dn^wn) ■ (4.13) 
consequently the classical elliptic theory combined with ()4.12p implies that 

/ \\Qn\\'^_i ,oa odt < K, uniformly with respect to (4-14) 

Jo 2+ ' 

From ([3ll) - ([3^ and ([iT2|) we also obtain that 

/ ||5tiUAr|P 3 , r,/, r,dt < K, uniformly with respect to A^. (4-15) 

Jo 2^ ' 

Step 2 (Passing to the limit N — t- oo) The central issues is how to take the limit in the 
nonlinear term D^fiiw^) x V x 0^,9(111^1^). 

From the Aubin-Lions compactness Lemma (the same arguments as in section 3) we can 
find a subsequence (wj^., Diy^^QWiy^,qiy^) and {w,z,q) such that when Nj — 00 we have: 



to AT - 


, * 

^ w 


weakly* in L°°{0,T; HI), 


(4.16) 


wn - 


^ w 


weakly in L^{id,T;Hl+^), 


(4.17) 


DnmWn - 


z 


weakly in L'^{0,T;HI), 


(4.18) 


dtWN - 


dtw 


weakly in L^{0,T; R-l^^'^), 


(4.19) 


QN - 


q 


weakly in L^{0,T; R-^+^^T^)), 


(4.20) 


Wn - 


w 


strongly in L^{0, T; H^) for ah s < 1 + 6*. 


(4.21) 



The goal is to prove 

DNfiWN X V X DNfiWN Agit) X V X kew strongly in L^(0, T; L^(T3)^), (4.22) 
Thus it remain to show that 

Dn,9Wn Agw strongly in L^{<d, T; L^). (4.23) 

In order to show ()4.23p . we compute directly the difference between Dn^qWj^i and Kqw as follows 

< IIA0 [wn - W) \\l2(o,T;L^) + \\DNfiW - ^ew\\L^{f),T;L^) 
< C\\wn - ■«^llL2(o,T;//2e) + WDn^OW - Agli; ||^2(o,j..£,2) 

(4.24) 
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Where we have used the second part of lemma BTTl Hence, using that ^ < 1, (j4.17p and the first 
part of lemma ()4.ip we deduce ()4.23p . Finally the result above, ()4.23p . combined with ()4.18p 
and the uniqueness of the limit, allows us to deduce 



V X Dn^wn -^"^ -X ^ew weakly in L2(0,r;L^). (4.25) 

Consequently we get ()4.22p and by using ()4.1ip and (|2.6p we obtain 



DNfiWN X V X DNfiWN ^8W X V X Aqw weakly in L'^{0,T; R-^^'^^ (Tsf). (4.26) 

These convergence results allow us to prove in the same way as in [5] that (if,g) is a dis- 
tributional solution to the mean rotational Navier-Stokes equations (|4.ip - (|4.4p . and (Agio, A^g) 
is a distributional solution to the rotational Navier-Stokes equations ()l.ip - ()1.2p . so we will not 
repeat it. 
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